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Abstract 

In this paper we present some recent results concerning the ex- 
istence, the stabihty and the dynamics of sohtons occurring in the 
nonhnear Schroedinger equation when the parameter h —?■ 0. 

We focus on the role played by the Energy and the Charge in 
the existence, the stability and the dynamics of solitons. Moreover, 
we show that, under suitable assumptions, the soliton approximately 
follows the dynamics of a point particle, namely, the motion of its 
barycenter q{t) satisfies the equation 

q{t) + VV{q{t)) = Hh{t) 

where 

sup\Hh{t)\^0 as h^O. 

Mathematics subject classification. 35Q55, 35Q51, 37K40, 37K45, 
47J35. 

Keywords: Soliton dynamics. Nonlinear Schroedinger Equation, or- 
bital stability, concentration phenomena, semiclassical limit. 

1 Introduction 

Roughly speaking a solitary wave is a solution of a field equation whose 
energy travels as a localized packet and which preserves this localization in 
time. 
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By soliton we mean an orbitally stable solitary wave so that it has a 
particle-hke behavior (for the definition of orbital stability we refer e.g. to 

Ref. [SlEllHllIilElEllES]). 

The aim of this paper is to review some recent results about the existence, 
the stability and the behavior of the solitary waves relative to the equation 

(1) 

ip (0, x) = ip{x) 

where (p{x) is a suitable initial data. 

In the first section we examine the case V = and h = 1 (see Ref. [3]). 
Under suitable assumption on W, there exists a stationary solution of the 
form tjj (t, x) = U{x)eT^'^^, where [/ is a radial function decaying at infinity 
which solves the equation 

- AU + W'{U) = 2ioU. (2) 

This solution is found by a constrained minimization method that involves 
two prime integrals of the motion: the Charge and the Energy. By a con- 
centration compactness argument it is proved that this stationary wave is 
stable, so this solution is a soliton. 

In the second part, we consider V ^ and h small (see Ref. ^). The 
stationary solution [/(x)e^'^* becomes 

^P{t,x) = u(^]e^ (3) 



where 

/S = 1 + f • (4) 

If /3 > 0, the stationary solution concentrate as /i — > 0. Also, we can give 
a precise estimate of the behavior of the Energy and the Charge. 

These estimates are the key ingredient to study the case V ^ 0. For h 
sufficiently small, a solution of ([I]) with initial datum ip{x) = [U (^^)] e^"^'^ 
is a soliton which travels like a point particle under the action of the potential 
V. In fact, in the last section of this review, we define the barycenter q{t) of 
a soliton (see Ref. as 



x\'ip{t, x)\'^dx 

q{t) = (5) 

\^{t, x)\'^dx 
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and we prove that it evolves approximatively like a point particle in a poten- 
tial V. More exactly, q{t) satisfies the Cauchy problem 



g(0) = go 

g(0) = V 
where 

sup \Hh{t)\ as /i 0. 

igR 

In the last years there are some result about existence and dynamics of 
soliton for the Nonlinear Schroedinger Equation (see, for example Ref. [6l 
El [ini [m 113 IISI ED [22]), in particular there are results which compare the 
motion of the soliton with the solution of the equation 

X{t) +VV{X{t)) = (6) 

for t G (0, T] for some constant T < oo. 

The result of Bronski and Jerrard [7] deals with a pure power nonlinearity 
and a bounded external potential. The authors have shown that if the initial 
data is close to f/(^^)e*"^ in a suitable sense then the solution iphit^x) of 



(dD satisfies for t € (0, T] 

■^\%ljh{t,x)\^ - j^^ \%IJh{t,x)\^dx!^ 5x(t) 



^ as /i ^ 0. (7) 



c 



Here 5x(t) is the Dirac "(5-function" , C^* is the dual of and X{t) satisfies 
the equation ([HD with X(0) = go, -^(0) = fo- 

In related papers of Keraani pTj [18] there are slight generalizations of the 
above result. Using a similar approach, Marco Squassina [22] and Alessan- 
dro Selvitella ^21j described the soliton dynamics in an external magnetic 
potential. 

Other results on this subject are in Ref. [Tni[II]. In Ref. (TU] the authors 
study the case of bounded external potential V . 

In Ref. [IT] the authors study the case of confining potential. They as- 
sume the existence of a stable ground state solution with a null space non 
degeneracy condition of the equation 

- Ar/^ + fifj^ + W'{r]^) = 0. (8) 

The authors define a parameter e which depends on fi and on other parame- 
ters of the problem. Under suitable assumptions they prove that there exists 
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T > such that, if the initial data tjj^{x) is very close to e*^°'*^^' "'''^~^^^°r]^g{x — 
ao) the solution x) of problem (Pi) with initial data ip^ is given by 

^{t, x) = e^p(*)-(^-«W)+*^(*),^^(,) (x - a{t)) + w{t) (9) 

with I Iz/i < e, p = — Vl^(a) + o(e^), a = 2p + 0(5^) with < t < j for e 
small. 

In our paper |1] we do not require the uniqueness of the ground state 
solution which is, in general, not easy to verify, and we formulate our result 
such that it holds for any time t. 

2 Main assumptions 

In all this paper we make the following assumptions: 

(i) the problem ([1]) has a unique solution 

G C°(M,ij2(R^)) nC^(R,L2(R^)) (10) 

(sufficient conditions can be found in Kato [16], Cazenave [9], Ginibre- 
Velo [13]). 

(ii) W : — > M is a function which satisfies: 

W{0) = W'{0) = W"{0) = (11) 

2N 

\W"{s)\ < cilsl"-^ + C2\s\f-HoT some 2<q<p<2* = — — -. (12) 

4 

W{s) > -c|s|^ c>0,2<z/<2 + — ands large (13) 

3so G M+ such that W{so) < (14) 

(iii) V : — > M is a function which satisfies the following assumptions: 

V{x) > 0; (15) 

\VV{x)\ < V{xf for |a;| > i?i > 1,6 G (0,1); (16) 
V{x) > for \x\ > Ri> l,a> 1. (17) 

(iv) the main assumption 

a > (18) 
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Let us discuss the set of our assumptions: 

The first assumption gives us the necessary regularity to define the barycen- 
ter and to prove that q(t) G C^(]R, M'^). The hypotheses on the nonhnearity 
are necessary in order to have a sohton type solutions. In particular, (fT2|) 
is a standard requirement to have a smooth energy functional, f|T4l) is the 
minimal requirement to have a focusing nonlinearity and (IT^ is necessary to 
have a good minimization problem to obtain the existence of a soliton. We 
require also that \^ is a confining potential (assumption (iii)). This is useful 
on the last part of this paper, to prove the existence of a dynamics for the 
barycenter. 

In our approach, the assumption a > is crucial. In fact, as we will 
see in Section 13.11 the energy of a soliton ip is composed by two parts: 
the internal energy Jh and the dynamical energy G. The internal energy is 
a kind of binding energy that prevents the soliton from splitting, while the 
dynamical energy is related to the motion and it is composed of potential 
and kinetic energy. We have that (see Section Hj) 

J,(V;) = /i^^-" 

and 

G(v^) = ll^lli. = /^^^ 

Then, we have that 

So the assumption a > implies that, for /i <C 1, G (ip) ^ (ip), namely 
the internal energy is bigger than the dynamical energy. This is the fact that 
guarantees the existence and the stability of the travelling soliton for any 
time. 
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2.1 Notations 

In the next we will use the following notations: 

Iie{z), Im(2;)are the real and the imaginary part of z 
B{xq, p) = {x E : |x — Xol < p} 
Sa = {u E : ||m||l2 = cr} 
r^ = {ueH': Mu) < c] 

m = m„2 := inf J{u) 

3 General features of NSE 



Equation ([T]) is the Euler-Lagrange equation relative to the Lagrangian den- 
sity 

1^ -m2 ...... ...^,^,..,2 (^g) 



where, in order to simplify the notation we have set 
Sometimes it is useful to write if) in polar form 



(20) 



(21) 



Thus the state of the system ip is uniquely defined by the couple of variables 
(m, S). Using these variables, the action S = J Cdxdt takes the form 



y I Vm|' + Wh{u) + Us + ^ \VS\' + Vix) ] 



dxdt 
(22) 



and equation ([1]) becomes: 

- ^Am + Wi{u) + {dtS + ^ \VS\^ + V{x)^ u = (23) 
dt [u^) + V ■ {u^VS) = (24) 
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3.1 The first integrals of NSE 

Noether's theorem states that any invariance under a one-parameter group 
of the Lagrangian imphes the existence of an integral of motion (see e.g. 
Gelfand-FominUg). 

Now we describe the first integrals which will be relevant for this paper, 
namely the energy and the "hylenic charge" . 

Energy The energy, by definition, is the quantity which is preserved by the 
time invariance of the Lagrangian; it has the following form 



\Vi^f + Wh{^) + V{x) 



dx. (25) 



EhW = 
Using fl?I]) we get: 

Ehiu,S) = J (^^\Vu\^ + Whiu)^dx + J Q|V5|' + V(x)^ u'^dx 

(26) 

Thus the energy has two components: the internal energy (which, 
sometimes, is also called binding energy) 

■Mu) = / (y I Vm|' + Wh{u)^ dx (27) 

and the dynamical energy 

Giu,S) = j (^^\VSf + V{x)^u^dx (28) 

which is composed by the kinetic energy | j iVSf u'^dx and the poten- 
tial energy j V{x)u'^dx. 

Hylenic charge Following Ref. [2] the hylenic charge, is defined as the 
quantity which is preserved by by the invariance of the Lagrangian 
with respect to the action 

For equation ([1]) the charge is nothing else but the norm, namely: 

niij) = [ \ijfdx = I u^dx. 



Momentum If ^ = the Lagrangian is also invariant by translation. In 
this case we have the conservation of the momentum 

Pjii;) = him J ^,^^dx, J = 1, 2, 3 (29) 
hence we have the first Newton law for the barycenter. 
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4 The case h = 1,V = 

In this section we present some results contained in Ref. We minimize 
the internal energy J{u) on the constraint {u G : ||tt||L2 = a} for some a 
fixed. If U is the minimizer and if 2a; is the Lagrange multiplier associated 
to U, ip{t,x) = f/(a;)e*'^* is a stationary solution of ([T]). 
We get the following result 

Lemma 1. Let W satisfy H^) . [T^) and [I4\ )- Then, 3 a such that W a > a 
there exists u G satisfying 

Jiu) = m„2 := inf JM, 

{v&m, |H|^2=a} 

with \ \u\\l2 = a. Then, there exist uj and u that solve with cu < and u 
positive radially symmetric. 

In order to have stronger results, we can replace f[T^ with the following 
hypothesis 

W{s) < < e < ^ for small s. (30) 

In this case we find the following results concerning the existence of the 
minimizer of J{u) for any cr. 

Corollary 2. // IfT^) . IfT^) and ^3^) hold, then for all a, there exists u G H^, 
with \ \u\\l2 = a, such that 

J{u)= inf J{v). 

{veH^ ,\\v\\ j^2=(t} 

In particular, for = 3 we have 

Corollary 3. Let N = 3. If ^^ and ^^ hold and W G , with W"'{fd) < 
0, then for all a, there exists u G with \ \u\\l2 = a such that 

J{u) = inf J{v). 

{veH'^ ,\\v\ \ j^2=(r} 

We sketch briefiy the steps of the proof for Lemma [H 
Step 1: If W satisfies ([HD then m^2:=inf5^ J{u) < 

Step 2: If W satisfies (I13p then m„2 > — oo, any minimizing Palais Smale se- 
quence Un is bounded in and the Lagrange multipliers Un associated 
to Un are bounded in M. 



8 



Step 3: Any minimizing Palais Smale sequence converges in to a minimizer. 

We point out that f|T^ is a fundamental requirement for the existence of a 
minimizer. In fact, if > 0, then by Pohozaev identity we can prove that 
f/ = is the unique radial solution of 

Concerning the stability of stationary solution we set 

S = {U{x)e'', e e S\ \\U\\l2 = a, J{U) = m,.} (31) 

Definition 4. S is orbitally stable if 

Ve, 35 > s.t. G iJ^(R^), inf || - u\h^ <6^ 

Vt inf II \'ijj{t, ■)\ - u\\h^ < e 
where ip is the solution of ([T]) with initial data ip. 

Using concentration compactness [HI [20] arguments we prove the follow- 
ing (see Ref. [5], Sect. 3) 

Theorem 5. Let W satisfy [W\) . [W\) and [T^ . Then S is orbitally stable. 

This variational approach can be successfully used to find stable solitary 
waves for the nonlinear Klein Gordon equation 

= W'm)^^. (32) 

Again, the crucial assumption to obtain solitons is f|T^ (see Ref. [1] for 
details) . 

We obtain a concentration result for a minimizer U crucial for this work 
(see Ref. [3]). 

Lemma 6. For any e > 0, there exists an R = R{e) and a 6 = 6{e) such 
that, for any u G J'^+^ ^ S^, we can find a point q = q{u) G such that 



1 

72 



u {x)dx < e. (33) 



We give a sketch of the proof. 

Proof. Firstly we prove that for any e > 0, there exists a 6 such that, for any 
u G we can find a point q = q{u) G and a radial ground state 

solution U such that 

\\u{x)-U{x-q)\\m<^- (34) 
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We argue by contradiction: if ([31]) do not hold, we can construct a minimizing 
sequence which not converge. At this point, given e, there exist a point 
q = q{u) G and a radial ground state solution U such that 

u{x) = U{x — q) + w and | |w| < Ce. (35) 

Now, we choose R such that 

^ / U^{x)dx < Ce (36) 

for all U radial ground state solutions. This is possible because if ?7 is a 
minimizer for J constrained on 5*0-, then there exists two constants C, R, not 
depending on U such that 

U{x) < Ce^l"'! for |x| >> R. 

By this fact we get the claim. □ 

We remark that, depending on the nonlinearity W, it is possible that the 
minimizer of the constrained problem is not unique. Anyway, by Lemma El 
R does not depend on the minimizer. 



5 The case h small enough 

We present now the main results contained in Ref. We recall some in- 
equalities which are useful in the following. Let it be 



We have 



and 



MI2 



Jbiu] 



2 + ^W{u)dx 
2 



^Af/3+2-2/3 



X 



/ X 



dx 



\V^v{i)\' + h''^--W{v{i))di 



h^P-- 1 i I (01' + W{v (0) d^ = /^^^-Vii 



(37) 
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We give now some results about the concentration property of the solu- 
tions ip{t, x) of the problem ([T]). Given K > 0, h > 0, we put 

ip{x) = i>{0,x) = Uh{0,x)ei^''^°'''^ 
with u,,iO,x)= [(U + w) (^)] 



B 



K 



(38) 



f/ is a minimizer of J constrained on Sa 

||;7 + ^i;||^2 = \ \U\\l2 = a and J{U + w) <m + Kh'' 

\\VSh{0,x)\\Lo^ < K for all h 

V{x)ul{0,x)dx < Kh^'^-^'^ 
Considering the set Bf^ as the admissible initial data set, we get 

Theorem 7. Assume V G and H^) . Fix K > 0, q eR^ . Let a> 0. 

For all e > 0, there exists R > and ho > such that, for any ipitjx) 
solution of (CP with initial data ipi^OjX) G Bj^ with h < Kq, and for any t, 
there exists qhif) G for which 

1 



\ip{t^ x)\ dx < e. 



(39) 



R^\B(gh(t),R/i/3) 

Here qh{t) depends on %p{t,x). 

We give the proof because it is simple and quite interesting. 

Proof. By the conservation law, the energy Ehiipitjx)) is constant with re- 
spect to t. Then we have 



Ehii;it,x)) 



Jh{uhiO,x)) + / ul{0,x) 



\VSh{0,x)\' 



+ V{x) 



dx 



< Mu,iO,x)) + ^a'h'^^ + Kh''^ 

where C is a suitable constant. Now, by rescaling, and using that '?/'(0,x) G 
5^''', we obtain 



Eh{^{t,x)) < h^^-'^JiU + w) + Ch^^ 

< h^l^~''{m + Kh'') + Ch^^ (40) 
= /i^^-"(m + ir/i" + Cr) = /i^''-"(m + /i"Ci) 
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where Ci is a suitable constant. Thus 

= Eh{ip{t,x))-G{^{t,x) 



\VShit,x)\ 



2 



+ V{x) 



Uh{t,x) dx 



< /i^''-"(m + /i"Ci) (41) 
because > 0. By rescahng the inequahty f l4T|) we get 

J{uh{t, h^x)) <m + h'^Ci (42) 

So, if a > 0, for h small by a simple argument and Lemma [6] we get the 
claim. □ 

Roughly speaking we have that Jhi'ip) — and G{ip) = h^^ and this 

is the key of the proof. 

To simplify in the following we take an initial data of the type 

^{x) = U (^-^\ e^-^ (43) 



where g, v are fixed. Obviously ^{x) E for some K. 

5.1 Existence and dynamics of barycenter 

We recall the definition of barycenter of 



x\il){t^ x)\^dx 

qnit) = ^ . (44) 

\tjj{t,x)\'^dx 



The barycenter is not well defined for all the functions ip G if (M ). Thus 
we need the following result: 

Theorem 8. Let ilj{t,x) be a global solution of (QP such that ilj[t,x) G 
C(R,H^{R^)) n (:7i(M,L2(R^)) with initial data V^(0,x) such that 



\x\\^p{0, x)\ dx < +00. 

Then the map qh{t) : R M^, given by ^ is C'^{R,R^) and it holds 
■ u\ ^^{hJ^N'4^{t,x)Vij{t,x)dx) 

*w = m:^. ■ 

V(x)V\i,{t,x)\'dx 

m:^. — ■ 
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We have the following corollary 

Corollary 9. Assume / fi^) and the assumptions of the previous theorem; 
then 



6 The final result 

6.1 Barycenter and concentration point 

We have two quantities which describe the properties of the travelling soliton: 
the concentration point q and the barycenter q. If we want to describe the 
particle-like behavior of the soliton the concentration point q seems to be 
the natural indicator: it localize at any time t the center of a ball which 
contains the larger part of the soliton. Unfortunately we do not have any 
control on the smoothness of q{t) (indeed q is nor uniquely defined). The 
barycenter, at the contrary, for a very large class of solutions has the required 
regularity, and the equation (j47j) is very similar to the equation of the motion 
we want to obtain. In this paragraph, we estimate the distance between the 
concentration point and the barycenter of a solution tp{t,x) for a potential 
satisfying hypothesis ffT^ and ffTTl) . say a con/inmg' potential. 

The assumption f[T7l) is necessary if we want to identify the position of 
the soliton with the barycenter. Let us see why. Consider a soliton ip{x) and 
a perturbation 



Even if ip {x) <ti ip{x), when d is very large, the "position" of ip{x) and the 
barycenter of ipd{x) are far from each other. In Lemma [T21 we shall prove 
that this situation cannot occur provided that (fTTI) hold. In a paper [5] in 
preparation, we give a more involved notion of barycenter of the soliton and 
we will be able to consider other situations. 

Hereafter, fixed > 0, we assume that ip{t,x) is a global solution of 
the Schroedinger equation (H]), ip{t,x) E C{R, H^) nC^{R, H'^), with initial 
data ip{0,x) G 5f with 5f given by fl38p . We start with some technical 
lemma. 

Lemma 10. There exists a constant L > such that 





(47) 



mt,x)\\i 



ipd(yX) = ip{x) + ip{x — d) , (i G M" 



N 



< 




V{x)ul{t,x)dx <L Vt G M, 
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The proof follows by estimating the energy. 
Lemma 11. There exists a constant Ki such that 

\qh{t)\ <Ki for tern. 

The proof follows by Lemma fTOl and by f|T7|) . Furthermore, we can choose 
i?2 such that 

u1{t, x)dx 

\x\>R2 ^ 1 

— 2 

x)dx 

Lemma 12. Given < e < 1/2, and R2 as in the previous lemma. 
We get 

1. sup^gjg < R2 + R{e)h^ < -R2 + 1; for all h <h and 5 <5 small 
enough. 

2. sup^gK \ qh{t) - qh{t)\ < aoa-i + 3-R3£ + R{e)h^, for any R3 > R2, and 
for all h small enough. 

The hardest part of the proof is the estimate of 

/RiVxB(0,i?3) - ^hit)) Ul{t, X)dx 



h 



/rjv '^\{t^ x)dx 



Using (|T7|) and the previous estimates we can conclude. 

We notice that Ri, R2 and R^ defined in this section do not depend on e. 

6.2 Equation of the travelling soliton 

We prove that the barycenter dynamics is approximatively that of a point 
particle moving under the effect of an external potential V{x). 

Theorem 13. Assume (i)-(iv). Given K > 0, let ij{t,x) e C(M,H^) D 
C^{M.,H^) be a global solution of equation ^j, with initial data in Bj^ , h < 
Hq. Then we have 

qhit)+VV{qh{t))=mXt) (48) 
with \ \Hh(t)\\ LOO goes to zero when h goes to zero. 
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Proof. We know by Theorem [HI that 



W{x)u1{t, x)dx 

m + '-^^ = 

u\{t^ x)dx 



(49) 



Hence we have to estimate the function 

/Riv[VV(x) - VV(g^(t))]M2(t,x)rfx 



Hhit) = [WV{qH{t)) - V\/(g^(t))] + 
By Lemma [TT] and Lemma Wl\ we get 



Irn t^hitj x)dx 



max 
i,j ^1,...,N 

\t\ <Ki+R2 + l 



9V(r) 



dxidxj 



(50) 



< M 



3L 



J + SRsE + R{e)h'' 



for any R3 > R2 and some M > 0. 
T'o GstiixictijC 

^^[Vr(x) - VV{qh{t))]ul{t,x)dx 
Irn 'f'ti^^ x)dx 
we spht the integral three parts. 

- VV{qnmnl{t,x)dx_ 



L3 



/rjv u\{t, x)dx 

/R^\B(gft(t),fl(£)fe/3) \'^V{x)\ul{t,x)dx ^ 

/jgjv "^hi^^ x)dx 



Irn u\{t^ x)d: 



X 



By the Theorem [7] and by Lemma fT2] we have L3 < Me. 
We have also 

Li<Ki + R2 + I. 



Using hypothesis (ITBi) we have 

L2 < Me + 



L 



.1-6 



(51) 



(52) 
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where b G (0, 1) is defined in flTB]) . Concluding we have 

b 

^i-f' + M{2 + 3R3)e + 2MR{e)h^. (53) 

At this point we can have sup^ \Hh{t)\ arbitrarily small choosing firstly 
sufficiently large, secondly e sufficiently small, and finally h small enough. □ 

Corollary 14. Let i){t,x) G C(M, iJ^) p C^{R^H^) he a global solution of 
equation (QP, with initial data ^p{x) = U{^^^)eTi^'^ where U is a radial min- 
imizer of J on Sa, Qq G M^, v G M^, and h < Hq. Then the bary center q 
satisfies the following Cauchy problem 

qh{t) + ^V{qh{t)) = Hh{t) 
qh{0) = qo 
QhiO) = V 

Proof. The initial data belongs to Bj^ for some K. We apply the previous 
results to obtain the equation for q. The initial data g(0) and g(0) are derived 
with a direct calculation. □ 



,,,,,, 3LM \ L 
\Hh{t)\ < 



7 The swarm interpretation 

In this section we present a different point of view on our problem. Although 
this approach is non rigorous, it provides some physical intuitions which 
are inspiring for a better understanding of the general framework. We will 
suppose that the soliton is composed by a swarm of particles which follow 
the laws of classical dynamics given by the Hamilton- Jacobi equation. This 
interpretation will permit us to give an heuristic proof of the main result. 
First of all let us write NSE with the usual physical constants m and h: 

= -^^^ + + ^(^)^- 

Here m has the dimension of mass and the Plank constant, has the di- 
mension of action. 

In this case equations (125]) and become: 

- I^Am + ]^Wl^u) + (^dtS + ^ \VS\^ + V{x)^ n = 0; (54) 
dt{u')+V-{u^^^ =0. (55) 
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The second equation allows us to interpret the matter field to be a fiuid 
composed by particles whose density is given by 



Pn = 



and which move in the velocity field 



(56) 
m 

So equation fISB]) becomes the continuity equation: 

dtpH + V ■ (p^v) = 0. 

If 

- + lw;Xu) « w, (57) 

equation ( 15^ can be approximated by the eikonal equation 

dtS + ^\VS\^ + V{x) = 0. (58) 
2m 

This is the Hamilton- Jacobi equation of a particle of mass m in a potential 
field V. 

If we do not assume equation ( 1551) needs to be replaced by 

dtS + — \VSf + V + Q{u) = (59) 
2m 

with 

_ - jnym) Au + ly^(^) 

The term Q{u) can be regarded as a field describing a sort of interaction 
between particles. 

Given a solution S{t, x) of the Hamilton- Jacobi equation, the motion of 
the particles is determined by Eq. (15Bi) . 

7.1 An heuristic proof 

In this section we present an heuristic proof of the main result. This proof 
is not at all rigorous, but it helps to understand the underlying Physics. 

If we interpret p-^ = u"^ as the density of particles then 

Pndx 
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is the total number of particles. By (jSH]), each of these particle moves as a 
classical particle of mass m and hence, we can apply to the laws of classical 
dynamics. In particular the center of mass defined in ([5]) takes the following 
form: 

^ J xmpndx ^ J xpndx ^ 
J mpT-idx J pndx 

The motion of the barycenter is not affected by the interaction between 
particles (namely by the term f l59|) ). but only by the external forces, namely 
by W. Thus the global external force acting on the swarm of particles is 
given by 

'f = - I VV{x)pndx. (61) 



Thus the motion of the center of mass q follows the Newton law 

f = Mq, (62) 

where M = J mp-^dx is the total mass of the swarm; thus by fl^Ul) . (1611) and 

), we get 

/ Wpndx _ J WVu^dx 
m J p-udx m J u^dx 

If we assume that the u{t,x) and hence pu(t,x) is concentrated in the 



m 

lat t 

point q{t), we have that 



VVu'^dx ^ VV{q{t)) j u^dx 
and so, we get 

mq{t) ^ -VV{q{t)) . 

Notice that the equation mq{t) = —W{q{t)) is the Newtonian form of 
the Hamilton- Jacobi equation fl58|l . 
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